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Abstract  Military  medical  evacuation  (MEDEVAC)  systems  respond  to  casualty  incidents  and  transport  the  most  urgent 
casualties  to  a  medical  treatment  facility  via  multiple  types  of  air  ambulance  assets.  Military  MEDEVAC  systems  are 
subject  to  an  uncertain  number  of  service  calls  and  each  service  call  demands  different  system  operations  depending  on 
type  and  the  priority  level.  Therefore,  military  medical  planners  need  an  air  MEDEVAC  asset  management  system  that 
determines  how  to  dispatch  multiple  types  of  air  assets  to  prioritized  service  calls  to  maintain  a  high  likelihood  of  survival  of 
the  most  urgent  casualties.  To  reach  this  goal,  we  propose  a  novel  binary  linear  programming  (BLP)  model  to  optimally 
locate  two  types  of  air  assets  and  construct  response  districts  using  a  dispatch  preference  list.  Additionally,  the  BLP 
model  balances  the  workload  among  assets  and  enforces  contiguity  in  the  first  assigned  locations  for  each  air  asset.  The 
objective  of  the  BLP  model  is  to  maximize  the  proportion  of  high-priority  casualties  responded  to  within  a  pre-specified 
time  threshold  while  meeting  performance  benchmarks  to  other  types  of  casualties.  A  spatial  queuing  approximation 
model  is  derived  to  provide  inputs  to  the  BLP  model,  which  thus  reflects  the  underlying  queuing  dynamics  of  the  system. 
We  illustrate  the  model  and  algorithms  with  a  computational  example  that  reflects  realistic  military  data. 

Keywords  Air  MEDEVAC  asset  management  •  casualty  triage  •  spatial  queuing  •  binary  linear  programming  •  military 
medicine 


1  Introduction 


Military  medical  evacuation  (MEDEVAC)  systems  save  lives  by  responding  to  casualty  incidents  and  transporting  the 
most  urgent  casualties  to  a  medical  treatment  facility  (MTF)  via  a  fleet  of  distinguishable,  geographically-dispersed  air 
ambulance  assets  (servers)  dedicated  to  the  response  of  casualties.  It  is  important  to  know  where  to  locate  distinct  classes 
of  air  assets  and  how  to  dispatch  them  in  order  to  maximize  the  likelihood  of  survival  of  the  most  urgent  casualties.  To 
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do  so,  it  becomes  necessary  to  triage  casualties  so  that  soldiers  with  the  most  severe,  life-threatening  injuries  -  where 
death  can  occur  if  response  is  not  immediate  -  receive  the  most  timely  treatment.  This  paper  examines  the  following 
hard  problems  found  in  military  MEDEVAC  systems:  1)  how  to  geographically  locate  two  classes  of  air  assets  and  2) 
how  to  construct  response  districts  for  all  air  assets  in  the  system.  The  problems  mentioned  above  are  hard  because  it  is 
not  obvious  where  to  locate  scarce  air  assets  or  when  to  dispatch  further  air  assets  (and  ration  closer  ones)  to  casualty 
incidents  in  response  time-dependent  military  medical  systems. 

Optimal  decision-making  within  MEDEVAC  systems  is  complicated  due  to  multiple  types  of  casualties  and  distin¬ 
guishable  air  assets.  The  ranking  of  evacuation  precedence  is  handled  through  a  straightforward  casualty  triage  scheme: 

—  CAT  A:  alpha  category  includes  urgent  and  urgent-surgical  casualties  that  need  to  be  treated  within  one  hour. 

—  CAT  B:  bravo  category  includes  priority  casualties  that  need  to  be  treated  within  four  hours. 

—  CAT  C:  charlie  category  includes  routine  casualties  that  need  to  be  treated  within  24  hours. 

The  fleet  includes  two  classes  of  air  assets  to  manage: 

—  U.S.  Army  helicopters  such  as  the  UH-60A/L  Black  Hawk  or  HH-60  MEDEVAC  platforms  are  normally  leveraged  to 
respond  to  casualties  based  on  their  speed,  range,  and  en  route  medical  capability. 

—  U.S.  Air  Force  helicopters  such  as  the  HH-60G  Pave  Hawk  are  also  available  to  respond  to  casualties  but  are  slower 
because  they  carry  heavy  weapon  systems  and  have  a  larger  flight  crew. 

When  effectively  employed  together,  both  Army  and  Air  Force  air  MEDEVAC  assets  can  provide  a  more  efficient  and 
responsive  military  MEDEVAC  system  for  overseas  combat  and  stability  operations. 

We  propose  a  novel  binary  linear  programming  (BLP)  model  that  optimally  locates  two  classes  of  air  assets  within 
a  military  medical  system  deployed  overseas  and  assigns  these  assets  to  casualty  locations  using  a  dispatch  preference 
list.  A  dispatch  preference  list  is  an  ordered  ranking  of  air  assets  for  each  casualty  location.  Additionally,  the  BLP  model 
balances  the  workload  among  assets  and  enforces  contiguity  amongst  the  first  assigned  location  for  each  air  asset.  The 
objective  of  the  BLP  model  is  to  maximize  the  proportion  of  high-priority  casualties  (CAT  A)  responded  to  within  a 
given  response  time  threshold.  The  BLP  model  incorporates  coverage  thresholds  for  low-priority  casualties  to  encourage 
prompt  service  for  all  casualties.  Accounting  for  the  queuing  dynamics  in  military  medical  systems  results  in  a  more 
realistic  optimization  model  in  which  air  assets  are  not  always  available  for  service.  Therefore,  we  also  propose  a  queuing 
model  to  derive  performance  statistics,  including  busy  probabilities  and  independence  correction  factors,  for  two  classes 
of  air  assets.  Realistic  busy  probabilities  and  removing  the  assumption  of  asset  independence  leads  to  more  effective 
resource  allocation  decision-making  for  military  medical  planners.  The  main  contribution  of  this  paper  is  a  BLP  model 
that  constructs  a  complete  set  of  MEDEVAC  response  districts  (modeled  as  dispatch  preference  lists)  of  air  assets  to 
casualties  that  explicitly  account  for  two  types  of  air  assets  and  multiple  casualty  types. 

This  paper  is  organized  as  follows.  Section  2  provides  a  literature  review  of  service  system  location  models  and  the 
probabilistic  nature  of  military  MEDEVAC  systems.  Section  3  introduces  the  approximate  Hypercube  algorithm  [16, 
14]  used  to  compute  queuing  system  performance  statistics.  Section  4  introduces  the  BLP  model  to  locate  air  assets 
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and  construct  dispatch  preference  lists.  We  present  computational  results  for  a  military  medical  system  in  section  5. 
Conclusions  and  future  work  are  presented  in  Section  6. 

2  Literature  Review 

A  number  of  existing  research  papers  focus  on  air  MEDEVAC  asset  optimization  models  for  military  medical  systems. 
Bastian  [2]  presents  a  multi-criteria  decision  analysis  (MCDA)  model  using  stochastic,  mixed-integer  goal  programming 
to  determine  the  minimum  number  of  air  assets  needed  at  each  MTF.  This  in  turn  maximizes  the  coverage  of  the  theater¬ 
wide  casualty  demand  and  the  probability  of  meeting  that  demand.  Simultaneously,  Bastian’s  model  also  minimizes  the 
maximal  MTF  evacuation  site  total  vulnerability  to  enemy  attack.  Fulton  et  al.  [13]  introduce  a  two-stage  stochastic 
optimization  modeling  framework  for  the  medical  evacuation  (ground  and  air)  of  casualties,  which  identifies  optimal 
casualty  evacuation  sites  and  MTF  sites  in  response  to  stochastic  demands  for  service.  While  Bastian’s  and  Fulton’s 
papers  also  account  for  the  stochastic  nature  of  military  medical  systems,  this  paper  is  distinct  because  it  leverages 
queuing  theory  to  derive  busy  probabilities  to  more  accurately  model  the  stochastic  nature  and  availability  of  air  assets. 
Bastian  and  Fulton  [5]  present  a  geospatial-based  decision-support  tool  that  identifies  which  air  asset  to  launch  in  response 
to  a  casualty  event  given  knowledge  of  terrain,  aircraft  location,  and  aircraft  capabilities.  Zeto  et  al.  [24]  examine  the 
pre-location  of  air  assets,  along  with  type  and  quantity,  to  maximize  the  theater-wide  coverage  while  balancing  air  asset 
reliability.  Bastian  and  Fulton  [5]  and  Zeto  et  al.  [24]  examine  location  and  dispatch  problems,  but  differ  from  the 
model  presented  in  this  paper  because  they  do  not  construct  a  dispatch  preference  list  and  complete  ordering  of  all  air 
assets  for  each  casualty  location.  Bastian  et  al.  [4,3]  examine  medical  evacuation  platforms  of  the  Future  of  Vertical  Lift 
(FVL)  program  in  support  of  brigade  operations.  The  specific  operations  analyzed  include  the  zero-risk  aircraft  ground 
speed;  different  aircraft  engine  trade-off  considerations;  and  the  effects  of  weaponizing  the  current  air  asset  fleet  will 
have  on  the  fleet’s  range,  coverage  radius,  and  response  time.  In  contrast  to  the  evaluation  of  future  platforms,  this 
paper  presents  a  model  that  can  be  used  at  the  tactical  or  operational  level  to  decide  where  to  locate  air  assets  and 
to  identify  how  to  dispatch  air  assets  to  casualty  incidents.  Bourna  [7]  develops  the  Medical  Evacuation  and  Treatment 
Capabilities  Optimization  Model  (METCOM)  that  considers  policy  effects  on  key  measures  of  effectiveness,  and  then 
optimizes  treatment  and  facility  capacities  for  given  casualty  flows.  This  paper  does  not  consider  the  impact  of  policy 
on  military  medical  systems  or  examine  the  impact  of  casualty  flows  on  facility  capacities.  Fulton  et  al.  [12]  present  a 
Monte  Carlo  simulation  to  evaluate  rules  of  allocation  and  planning  considerations  for  Army  air  ambulance  companies 
during  major  combat  operations  (MCO).  Rather  than  incorporating  a  simulation  model,  this  paper  uses  queuing  theory 
and  binary  linear  programming  to  model  and  examine  military  medical  systems. 

There  has  been  much  research  on  civilian  emergency  medical  service  (EMS)  systems  which  share  similar  characteristics 
as  military  medical  systems.  There  are  many  ambulance  location  models  that  focus  on  “covering”  patient  demand  based 
on  response  time  thresholds  that  are  (usually)  nine  minutes  from  dispatch.  Many  of  the  models  are  extensions  and 
variations  on  the  seminal  paper  by  Church  and  ReVelle  [10]  that  introduces  the  Maximal  Covering  Location  Problem 
(MCLP)  as  a  basic  facility  location  model  for  EMS  systems.  An  extension  to  MCLP  is  the  Maximal  Expected  Covering 
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Location  Problem  (MEXCLP)  developed  by  Daskin  [11],  which  attempts  to  maximize  the  expected  number  of  calls 
covered  in  a  given  amount  of  time,  and  makes  the  following  assumptions:  servers  operate  independently,  each  server 
has  the  same  busy  probability,  and  busy  probabilities  do  not  depend  on  server  location.  The  model  in  this  paper  lifts 
Daskin’s  last  two  assumptions.  ”  Assets  of  different  types  may  be  used  in  ways  that  result  in  large  differences  in  the  busy 
probabilities  between  the  asset  types.  Likewise,  there  are  often  “hot-beds”  of  activity  within  military  systems  (e.g.,  near 
an  enemy  stronghold)  that  would  naturally  result  in  some  assets  to  have  larger  busy  probabilities  due  to  their  frequent 
selection  as  the  responding  assets.  This  paper  addresses  both  of  these  issues.”  Batta  et  al.  [6]  further  extend  MEXCLP 
by  introducing  the  Adjusted  MEXCLP  (AMEXCLP)  to  lift  the  three  previously  mentioned  assumptions  of  MEXCLP 
by  use  of  Hypercube  correction  factors  (see  below)  for  one  type  of  server  and  a  single  casualty  type.  Our  model  extends 
the  work  by  Batta  et  al.  [6]  by  examining  more  than  one  type  of  server,  multiple  casualty  types,  and  probabilistic  travel 
times.  The  papers  thus  far  do  not  consider  distinguishable  types  of  ambulances  as  we  do  in  this  paper. 

ReVelle  and  Marianov  [23]  and  Marianov  and  ReVelle  [19]  examine  the  location  of  multiple  types  of  servers — fire 
engines  and  fire  trucks — to  maximize  the  number  of  calls  requiring  both  fire  engines  and  fire  trucks.  Revelle  and  Marianov’s 
work  assumes  that  one  server  of  both  types  is  required  for  each  call  for  service,  where  this  model  dispatches  one  server  to 
each  call  for  service.  Mandell  [18]  also  considers  multiple  types  of  medical  units  but  does  not  consider  multiple  casualty 
types.  Marianov  and  Serra  [20]  introduce  a  model  that  maximizes  the  demand  covered  when  the  customer  does  not  have 
to  wait  in  line  (due  to  congestion)  with  more  than  a  prespecified  number  of  other  customers.  In  contrast  to  Marianov 
and  Serra’s  work,  there  is  a  zero  length  queue  in  this  model.  McLay  [21]  introduces  a  model  that  locates  two  types  of 
ambulances  while  responding  to  multiple  types  of  customers  but  does  not  assign  response  zones  to  ambulances.  This 
paper  extends  the  work  of  McLay  [21]  to  construct  dispatch  preference  lists  modeled  as  contiguous  response  zones  while 
also  balancing  asset  workloads. 

Spatial  queuing  models  have  been  widely  used  to  describe  the  underlying  dynamics  of  military  medical  and  civilian 
EMS  systems.  The  exact  and  approximate  Hypercube  models  by  Larson  [15, 16]  are  the  most  well-known  spatial  queuing 
models.  The  Hypercube  model  assumes  a  multi-server  queuing  system  with  indistinguishable  servers.  Jarvis  [14]  extends 
the  Hypercube  queuing  model  with  an  approximation  algorithm  to  incorporate  dependencies  on  casualty  call  type  for 
service  times  with  servers  that  are  distinguishable  by  location.  Chelst  and  Barlach  [9],  Larson  and  McKnew  [17],  and 
Budge  et  al.  [8]  consider  other  extensions  to  the  Hypercube  model.  Several  papers  use  Hypercube  model  outputs  as  BLP 
model  inputs  (see  Batta  et  al.  [6],  McLay  [21],  and  Ansari  et  al  [1]). 

This  work  adds  to  the  literature  by  introducing  a  new  BLP  model  that  optimally  locates  two  classes  of  air  assets 
and  assigns  assets  to  casualty  locations  in  a  preference  list.  This  work  also  includes  an  extension  to  the  approximate 
Hypercube  queuing  model  that  accounts  for  multiple  casualty  types  to  derive  system  statistics  that  are  used  as  input 
to  the  BLP.  The  new  BLP  model  in  this  work  better  captures  the  military  medical  evacuation  problem  studied  than 
previous  models.  For  example,  recently  military  medical  planners  and  senior  decision-makers  must  manage  two  distinct 
classes  of  air  assets  (both  U.S.  Army  and  U.S.  Air  Force)  to  respond  to  triaged  casualty  events.  Further,  contiguous 
response  zones  are  a  realistic  detail  of  this  model  that  reflect  the  practice  of  a  determining  dedicated,  compact  area  of 
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responsibility  for  helicopter  pilots.  Accounting  for  the  queuing  dynamics  in  military  medical  systems  results  in  a  more 
realistic  optimization  model,  in  which  air  assets  are  not  always  available  for  service. 


3  Queuing  Model 

In  this  section,  we  introduce  a  queuing  model  that  computes  system  performance  statistics  to  reflect  the  complex  dynamics 
of  military  medical  systems.  Air  assets  are  not  always  available  for  service  and,  therefore,  computing  busy  probabilities 
allows  for  more  realistic  air  asset  management  in  military  systems.  Servers  do  not  act  independently  of  one  another  in 
service  systems;  thus,  the  queuing  model  in  this  section  also  develops  independence  correction  factors  to  account  for  the 
assumption  of  server  independence.  A  list  of  the  symbols  used  in  this  paper  is  given  next,  followed  by  a  discussion  of  the 
system  dynamics  and  assumptions  of  the  queuing  model. 

sa  (sb)  ~  Number  of  type  A  (type  B)  air  assets  in  the  system  with  sa  +  sb  =  s, 

I A  (Ib)  ~  Set  of  all  potential  type  A  (type  B)  MTF  locations, 

J  -  Set  of  all  casualty  locations, 

C  -  Set  of  casualty  types  that  can  arrive  partitioned  into  four  classes  Ca,  Cab,  Cba,  Cb, 

Ca  (Cb)  C  C  -  Set  of  casualties  that  require  type  A  (type  B)  air  asset  evacuation, 

Cab  (Cba)  CO-  Set  of  casualties  that  prefer  type  A  air  asset  evacuation  over  type  B  air  evacuation  (prefer  type  B 
over  type  A), 

Ac  -  Call  arrival  rate  of  casualty  type  c,  for  c  £  C,  with  Ac  =  J  Afo>  anc^  Afo  the  ca^  arrival  rate  of  a  casualty  at 
location  j  of  type  c,  for  j  £  J,  c  £  C , 

A  -  System- wide  total  call  arrival  rate  with  A  =  ceC  Afo! 

T-jc  (TijC)  ~  Service  time  when  type  A  (type  B)  air  asset  location  i  responds  to  casualty  of  type  c  at  location  j ,  for 

c  £  C ,  j  £  J ,  i  £  Ia(Ib) , 

ta  (tb)  -  System-wide  average  type  A  (type  B)  air  asset  service  time, 

9C  -  Minimum  performance  benchmark  for  the  fraction  of  casualty  type  c  calls  that  must  be  responded  to  in  a  prespecified 
time  threshold,  for  c  £  C, 

PA  (pB)  ~  Traffic  intensity  for  type  A  (type  B)  air  assets, 

rf  (rf)  -  Type  A  (Type  B)  air  asset  location  i  busy  probability,  for  i  £  Ia(Ib), 
rA  (rB)  -  System- wide  average  type  A  (type  B)  air  asset  busy  probability, 

cijck  “  Ordered  dispatch  preference  list  where  cijCk  is  the  kth  preferred  type  A  air  asset  for  casualty  of  type  c  at  location 
j ,  for  j  £  J,  c£  Ca  U  Cab,  k  =  1, . . . ,  sa  or  for  c  £  Ca  U  Cba,  k  =  sB  +  1,  •  •  • ,  sB  +  sa, 

■fijck  (fijck)  '  The  dispatch  probability  of  type  A  air  asset  location  Uni  a  (type  B  air  asset  i  £  Ib)  to  casualty  j  £  J  oi 
type  c  £  C  as  the  kth  =  1, . . . ,  sa  (sb)  priority, 

PA  (PB)  -  Loss  probability  that  k  type  A  (type  B)  air  assets  are  busy, 
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R.fjc  ( RBC  )  -  Fraction  of  time  type  A  (type  B)  air  asset  location  i  can  respond  to  casualties  at  location  j  of  type  c  in  a 
prespecified  time  threshold,  for  i  G  Ia{Ib ),  j  G  J,  c  €  C, 

Qa(sa,  pA ,  k)  or  Qb  (sb,  pB  ,k)-  Hypercube  correction  factor  as  a  function  of  the  number  of  air  assets  sa  or  sb,  traffic 
intensity  pA  or  pB ,  and  priority  k,  for  k  =  0, sa  —  1  or  sb  —  1, 

Nj  -  A  neighbor  incident  matrix  used  for  contiguity,  i.e. ,  the  set  of  casualty  locations  that  are  adjacent  to  casualty 
location  j  in  the  usual  sense,  for  j  G  J. 


Our  queuing  model  is  an  extension  to  the  approximate  Hypercube  queuing  model  given  in  Jarvis  [14]  that  preserves  the 
location  of  assets.  The  extension  considers  two  types  of  air  assets — type  A  and  type  B — and  multiple  casualty  types  in 
a  military  MEDEVAC  system.  Consider  a  system  with  sa  type  A  air  assets,  sb  type  B  air  assets,  J  casualty  locations, 
and  C  casualty  types.  Classification  by  triage  leads  to  four  mutually  exclusive  casualty  types:  (1)  casualties  that  require 
type  A  air  asset  response  denoted  Ca,  (2)  casualties  that  prefer  type  A  air  asset  response  over  type  B  air  asset  response 
denoted  Cab ,  (3)  casualties  that  prefer  type  B  air  asset  response  over  type  A  air  asset  response  denoted  Cba,  and  (4) 
casualties  that  require  type  B  air  asset  response  denoted  Cb ■  Note  that  C  =  CaCCabCCbaCCb-  Casualties  at  location 
j  G  J  of  type  cGC  arrive  according  to  a  Poisson  process  with  rate  A  jc  independent  of  other  casualty  locations  or  types. 
The  total  arrival  rate  of  casualties  in  the  system  is  A  =  Y2je.j  SceC  Ajc-  Exactly  one  air  asset  is  assigned  to  each  casualty 
for  evacuation  unless  all  air  assets  are  busy,  in  which  case  a  non-traditional  evacuation  air  asset  will  be  called  upon  for 
assistance.  In  military  medical  systems  there  is  a  zero  length  queue  for  service.  Note  that  in  this  section  we  develop 
formulas  for  the  type  A  air  assets.  The  equations  for  the  type  B  air  assets  are  omitted  but  computed  analogously. 


In  military  medical  systems,  air  assets  do  not  act  independently  from  one  another  and  thus  an  optimization  model  to 
locate  and  dispatch  air  assets  should  account  for  asset  dependencies.  As  in  Jarvis  [14],  the  type  A  air  asset  independence 
correction  factors  Qa(sa 5  PA,j )  quantify  the  correction  to  the  probability  of  obtaining  j  busy  type  A  air  assets  followed 
by  an  available  type  A  air  asset  when  assuming  that  air  assets  operate  independently.  Note  that  Paa  is  the  probability 
that  all  sa  type  A  air  assets  are  busy. 


nA(=A  A,.,  forfc_01  i 

C  [SA,P  /  ,  /  ■_  i,\\ _  t/i  r>A  \k-r-t  _4/i  r>4\\’  tor  ^  0,1,...,  Sa  1- 


j=k 


(j  -  A:)!sa!(1  -  PA)k(  1  -  PA(  1  -  PA))  '■ 


(1) 


It  is  desirable  to  know  the  probability  of  assigning  specific  air  assets  to  specific  customers.  Denote  the  dispatch 
probability  of  type  A  air  asset  location  Uni  a  to  casualty  j  G  J  of  type  cGC  as  the  kth  =  1, . . . ,  sa  priority  as  ,fAck- 
For  type  A  air  assets,  fAck  is  computed  differently  for  casualty  types  Ca  and  Cab  than  casualty  type  Cba  since  type 
A  air  assets  can  only  be  dispatched  to  casualties  of  type  Cba  if  all  type  B  air  assets  are  busy  in  service.  The  remaining 
system  statistics  depend  on  fAck ■  Assume  there  is  a  known  dispatch  preference  list,  and  let  a,jck  equal  the  air  asset  that 
is  kth  preferred  to  respond  to  a  casualty  of  type  c  at  location  j.  For  example  if  a\A2  =  3,  then  air  asset  3  is  the  second 
preferred  air  asset  to  respond  to  a  casualty  of  type  A  at  location  1.  In  the  computation  of  fAck,  rA  is  the  busy  probability 
for  the  type  A  air  asset  location  UiiIa  ■  Note  that  fAck  =  0  for  c  G  Cjj  since  type  B  air  assets  are  required  to  respond  to 
casualties  of  type  Cb- 
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k- 1 

.fijck  ~  QA(sa,  PA,  k  -  1)(1  -  rf)  rfjcl ,  for  c  £  CA  U  Cab,  j  £  J,  k  =  1, . . . ,  sA,  i  £  I  a  (2) 

Z=1 

fe-1 

/iicft  ~pFBQA(sA’PAik-  1)(1  -  rf)  r£jcl,  force  CBA,  j  &  j,  k  =  sB +  1,.  ..,sB  +  sA,  i  £  Ia  (3) 

1=sb  T1 

The  dispatch  probabilities  computed  in  equations  (2)  and  (3)  depend  on  three  values:  1)  correction  of  the  assumption 
of  independence  among  air  assets  Qa(sa,  pA ,k  —  1),  2)  the  probability  that  the  specific  type  A  air  asset  is  available 
(1  —  rf),  and  3)  the  probability  that  the  k  —  1  more  preferred  air  assets  are  busy  Yl^i  rajcl-  Note  that  a,jci  =  i  is  a 
dispatch  preference  list,  so  the  formula  for  fAck  is  multiplying  together  the  busy  probabilities  of  the  k—  1  more  preferred 
air  assets.  Lastly,  equation  (3)  includes  the  PfB  term  because  a  type  A  air  asset  is  only  assigned  to  a  Cba  type  casualty 
once  all  sB  type  B  air  assets  are  busy. 


The  computation  of  the  type  A  air  asset  busy  probabilities  rf  is  given  by  equation  (4). 


rf  = 


V/ 


(1  +  ^A) 


,  i  £  I  A 


(4) 


In  equation  (4),  VA  is  the  rate  at  which  type  A  air  asset  i  is  assigned  to  a  call  for  service  and  is  computed  by  equation 
(5): 


SA 

ea  =  E 


E 

ceC AUC AB  j  : 


E 


i(sA,pA,  k  -  1)  P  rf. 


+  p* 


B+SA 

E 

=  »B  +  1 


E  E 

CBA  3  :  ajck  z 


xjorfcQA{sA,pA,k-l) 


(5) 


for  i  £  I  a-  Equation  (5)  is  comprised  of  two  pieces  separated  by  addition.  The  first  piece  captures  the  rate  at  which  type 
A  air  asset  i  is  assigned  to  casualty  types  Ca  and  Cab-  Type  A  air  asset  i  is  assigned  to  casualties  only  when  the  more 
preferred  air  assets  are  busy,  captured  by  JljLi1  rfjcl  ■  The  first  summation  shows  that  type  A  air  asset  i  can  be  any  of 
the  first  preferred  air  assets.  The  third  summation  is  over  the  casualty  location  j  only  if  the  dispatch  preference 

list  maps  type  A  air  asset  i  to  casualty  location  j  as  the  kth  preferred  asset.  The  second  piece  of  equation  (5)  captures 
the  rate  at  which  type  A  air  asset  i  is  assigned  to  casualty  types  Cba-  Note  that  the  preference  index  k  start  at  sB  +  1 
because  type  A  air  assets  cannot  be  assiged  to  the  first  sB  spots  of  the  dispatch  list  for  casualty  type  Cba-  Also,  PfB  is 
included  in  this  second  piece  because  a  type  A  air  asset  is  only  assigned  to  a  Cba  type  casualty  once  all  sB  type  B  air 
assets  are  busy. 


The  average  service  time  for  type  A  air  assets  ta  is  given  by  equation  6. 

tA  =  TijJijck/( 1  -  pAA)  (6) 

j(zJ  cGC  i= 1  k— 1 

In  equation  (6) ,  The  average  service  time  is  a  product  of  the  service  times  rfc  and  the  normalized  dispatch  probabilities 
fijck/ (1  —  PAa)-  The  average  service  time  is  weighted  by  the  call  arrival  rate  to  casualties  (AJC/A). 


Benjamin  C.  Grannan,  Nathaniel  D.  Bastian,  and  Laura  A.  McLay 


The  approximate  algorithm  is  run  until  the  maximum  change  in  air  asset  busy  probabilities  falls  below  a  pre-specified 
e.  The  approximate  algorithm  contains  an  initialization  step  where  the  closest  air  asset  is  assumed  always  available. 

Given: 

1.  Call  arrival  rates  A jC,  for  j  G  J,c  G  C 

2.  Type  A  air  asset  service  times  tac,  for  i  €  Ia,3  G  J,c  £  Ca  U  Cab  U  Cba 

Initialize: 


rf  =  Yl  Xjct£c,  i  G  iA  (7) 

c  j:i  is  the  closest 

tA  =  X/ ( Xi*/X)Tijc  (8) 

j£j  cGC  k= 1 

Iterate: 

Step  1.  Compute  <5a(sa,  pA,  k)  for  k  =  0, 1, _ ,  —  1  using  equation  (1).  Use  pA  —  A ta/sa- 

Step  2.  For  i  =  1, . . . ,  sa,  compute  rf  with  equation  (4),  where  VA  is  given  by  equation  (5). 

Step  3.  Stop  if  maximum  change  in  rA  is  less  than  a  tolerance. 

Step  4-  Else,  compute  ta  by  equation  (6),  where  /^cfe  is  given  by  equations  (2)  and  (3). 

Step  5.  Return  to  Step  1. 

To  compute  the  type  B  correction  factors  and  average  air  asset  busy  probabilities  also  used  as  inputs  to  the  BLP 
model,  the  analogous  type  B  approximate  algorithm  is  run  with  the  corresponding  type  B  air  asset  equations.  The 
following  differences  are  necessary  to  derive  the  type  B  equations.  Type  B  air  assets  are  preferred  first  when  responding 
to  casualties  of  type  c  G  Cb  U  Cba  as  priority  k  =  1 , ,sB-  The  type  B  air  assets  responding  to  casualties  of  type 
c  G  Cab  are  preferred  after  type  A  air  assets  in  the  dispatch  list  as  priority  k  =  sa  +  1, . . . ,  sa  +  sB-  The  type  A  air 
asset  independence  correction  factors  Qa(sa,  pA,k)  and  average  air  asset  busy  probabilities  rA  =  Y^ieiA  rtlsA  are  used 
as  input  to  the  BLP  model  formulated  in  the  next  section.  Likewise,  the  type  B  air  asset  independence  correction  factors 
QB (sb,  pB ,k)  and  average  air  asset  busy  probabilities  rB  =  Y^ieiB  r? /sb  are  also  used  as  input  to  the  BLP  model.  The 
next  section  introduces  an  optimization  model  that  utilizes  the  queuing  correction  factors  and  dispatch  probabilities. 

The  outputs  of  the  spatial  queuing  approximation  can  be  used  to  compute  several  inputs  to  the  BLP  model  derived 
in  the  next  section.  The  coefficients  hAck  and  hBck  found  in  the  objective  function  and  constraints  of  the  BLP  in  the 
next  section  represent  the  fraction  of  casualties  of  type  c  at  location  j  that  air  asset  i  can  reach  in  a  prespecified  time 
threshold  as  the  fctli  preferred  asset  to  send,  weighed  by  the  demand  at  location  j  (captured  by  the  Ajc/Ac  term).  These 
coefficients  utilize  dispatch  probabilities  (see  (2  and  3)  when  assuming  that  busy  probabilities  of  assets  of  the  same  type 
are  equal.  The  coefficients  for  type  A  assets  hf-ck  are  as  follows: 


htjck  =  (xjc/Xc)R?jC  QA{sA,pA,k)(l  -  rA)  (rA)k  1  ,i  G  IA,j  G  J,c  G  CA  UCAb,  k  =  1,...,sa 


(9) 


hfjck  =  ixjc/xc)  RfjcPFB  QA{.SA,pA,k)(l-rA){rA)k  1  ,  i  G  I  a,  j  G  J,  c  G  CBa,  k  =  sB  +  1, ...,  sB  +  sa- 


(10) 
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(ii) 


These  coefficients  take  two  sources  of  uncertainty  into  account:  (1)  uncertain  asset  availability  reflected  in  the  dispatch 
probabilities  and  correction  factors  and  (2)  uncertain  travel  times  that  lead  to  imperfect  coverage,  (reflected  by  the  terms 
Rfjc  that  represent  the  fraction  of  type  c  calls  that  can  be  reached  in  the  fixed  response  time  threshold). 

In  a  similar  manner,  the  coefficients  gAjck  and  gfck  represent  the  asset  workloads  for  type  A  and  type  B  air  assets, 
respectively,  and  the  coefficients  for  type  A  assets  gAjck  are  defined  in  the  following  equations: 


dfjck  =  (Atc/Ac)rA  QA(sA,pA,k)(l  -  rA)  (rA)k  1  ,  i  £  IA,  j  &  J,  c  £  CA  U  Cab,  k  =  1, ...,  sA 


gfjck  =  (Aic/Ac)  rAcPfB  Qa(sa,  pA,  k)(  1  -  rA)  ( rA ) 


i  i  G  Iai  j  €  J,  c  £  CBA ,  k  —  sB  +  1,  +  sB 


(12) 

(13) 

(14) 


where  tA,  and  t$c  are  the  service  times  associated  with  type  A  and  type  B  assets,  respectively.  These  coefficients  are 
used  to  balance  the  workload  among  the  servers  of  the  same  type,  and  therefore,  they  also  reflect  the  queuing  dynamics. 


4  Biliary  Linear  Programming  Model 

We  next  present  the  BLP  model  used  to  locate  two  classes  of  air  MEDEVAC  assets  and  construct  dispatch  preference 
lists  while  balancing  the  workload  among  assets  and  enforcing  contiguity  amongst  the  first  assigned  location  for  each  air 
asset.  The  objective  of  the  BLP  model  is  to  maximize  the  proportion  of  high-priority  casualties  (CA  calls)  responded  to 
within  a  pre-specified  time  threshold.  The  BLP  model  also  incorporates  coverage  thresholds  for  low-priority  casualties  to 
provide  prompt  service  for  all  casualties. 

To  describe  the  BLP  model  we  first  introduce  the  decision  variables.  The  first  two  sets  of  binary  variables  yf  and  yf 
indicate  the  set  of  MTF  locations  utilized:  yf  =  1  if  type  A  MTF  location  i  £  IA  is  utilized  and  0  otherwise.  Likewise, 
yf  =  1  if  type  B  MTF  location  i  £  IB  is  utilized  and  0  otherwise. 

The  next  sets  of  binary  variables  xfck  and  xfck  capture  the  dispatch  preference  fist  assignment  of  type  A  and  type 
B  air  assets  to  casualties.  These  decision  variables  are  defined  to  reflect  the  response  protocols  to  casualties  of  different 
types.  For  example,  only  type  A  air  assets  can  be  dispatched  to  casualties  of  type  CA,  and  type  A  air  assets  must  be 
dispatched  to  casualties  of  type  CAB  before  type  B  air  assets  are  dispatched  as  backup  coverage  to  casualties  of  type  CAB . 
Likewise,  only  type  B  air  assets  can  be  dispatched  to  casualties  of  type  CB,  and  type  B  air  assets  must  be  dispatched  to 
casualties  of  type  CBA  before  type  A  air  assets  are  dispatched  in  secondary  support  to  casualties  of  type  CAB .  Therefore, 
xtjck  =  1  if  type  A  air  asset  i  £  IA  is  the  kth,  k  =  1, . . .  ,sA,  preferred  air  asset  to  respond  to  a  casualty  at  location 
j  £  J  of  type  c  £  CA  U  CAB ,  or  if  type  A  air  asset  location  iinIA  is  the  kth,  k  =  sB  +  1, . . . ,  sB  +  sA,  preferred  air  asset 
to  respond  to  a  casualty  of  type  c  £  CBA ,  and  0  otherwise.  Likewise,  xfck  =  1  if  type  B  air  asset  location  iinIB  is  the 
kth,  k  =  1, . . . ,  sB,  preferred  air  asset  to  respond  to  a  casualty  at  location  j  £  J  of  type  c  £  CB  U  CBA ,  or  if  type  B  air 
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asset  location  Unis  is  the  kth,  k  =  SA  +  l,---, sa+sb,  preferred  air  asset  to  respond  to  a  casualty  of  type  c  G  Cab,  and 
0  otherwise. 

Now,  we  formally  state  the  binary  linear  program. 


max 

cECa 

EEE  htjckxtjck 
j€J  ^G/a  k— 1 

(15) 

subject  to 

El  1 

*£/a 

v  j  G  J, 

c  G  Ca  U  Cab,  k  =  1, . . 

• ,  sa 

(16) 

^  1  ^i/jck  i 

»£/a 

v  j  G  J, 

c  G  Cba,  k  =  sb  +  1, .  • 

• ,  sb  +  sa 

(17) 

E  = 1 

^G/b 

Vj  G  J, 

c  G  Cb  U  Cba,  k  =  1, . . 

• ,  sb 

(18) 

El  ^V/cfc _  i 

ieiB 

SA 

Vj  G  J, 

c  G  Cab,  fc  =  sa  +  1,  •  • 

. ,  sa  +  sb 

(19) 

EE  ^ijck^ijck  >  ec  Vc  e  cba 

Eg/b  k— 1  /  / 


E  xt>ck  —  vt  Vj  G  </,  c  G  C  a  U  Cab,  i  G  I  a 

fc=i 
SB+SA 

E  xtick  =  y?  Vj  G  J,  C  G  Cba,  i  G  Ct 

fc  =  SB+l 

SB 

E  =  2/f  Vj  G  J,  c  G  Cb  U  Cba  *  G  /b 

fc  =  l 
SA+SB 

E  xijck  =  y ?  Vj  G  J,  c  G  CAb,  i  &  Ib 

k—SA~\~  1 

«  SA  \  /  SA+SB  \  \ 

EE  h'ijck'X'ijck  ME  E  a  ^ijck^ijck  ))>0c  Vc  G  Cab 

i^lIa  k= 1  /  V^G-Zb  &=sa  +  1  /  / 

«SB+SA  \ 

E  E  P^Bhtjckxtjck  ) 

£G/a  fc=ss  +  l  / 

SB 

EEE  ^'ijck'^ijck  >9C  VcgCb 

j€J  iG-Zs  fe=l 

E 

»e/A 

E  v? =  ss 

iG/n 

Syc*  fs  yt  Vi  G  /a,  j  G  J,  c  G  Ca  U  Cab,  k  =  1, . . . ,  sa 
xtjck  —  yf  Vi  G  /a  ,  j  G  J,  c  G  Cba  ,  /c  =  sb  +  1  ,  •  •  •  >  sb  +  sa 
<  yf  Vi  G  Ib ,  j  €  J,  c  G  Cb  U  Cba,  k  =  1, . . . ,  sb 
xfjck  —  yf  Vi  G  /bj  j  G  J,  c  G  Cab,  fc  =  sa  +  5  sa  +  sb 

SA  SB+SA 

EEE  Sijck^ijck  EE  E  PsB9tjckxtjck  <  (rA  +  SA)yf  Vi  G  /A 

jG  J  cECba  Zc=sb-I-1 

SA+SB 

EEE  dijck^ijck  EEE  Ptafjckxfjck  <  {rB  +  8B)yf  Vi  els 

jEJ cECbUCba  fc= 1  je+ cSCab  fc=SA  +  l 


E 

E 

E 

jEJ  cECaUCab 

k= 1 

E 

E 

sb 

E 

(20) 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 
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'y  '  xtj'c i  —  xtjci  V*  G  Ia,c  £  Ca  u  Cab ,  j  £  J 

feNj 

(35) 

^2  xij'c  1  -  xijc i  v*  e  Ib,  C  £  CB  u  Cba,  j  G  J 

j’ZNj 

(36) 

xtjck  G  {0, 1}  Vi  £  I a,  j  £  J,c  £  C,k  =  1, . . . ,  s a  +  sb 

(37) 

xijck  G  {o,  1}  Vi  £  IB,j  £  J,c  £  C,k  =  1, sa  +  Sb 

(38) 

yf  £  {0, 1}  Vi  £  I A 

(39) 

yf  £  {0, 1}  Vi  e  IB 

(40) 

The  objective  function  in  (15)  maximizes  the  expected  coverage  of  Ca  casualties.  In  the  BLP  model,  only  type  A  air 
assets  respond  to  Ca  casualties.  Cab  casualties  prefer  type  A  air  assets  due  to  operational  characteristics  such  as  in¬ 
flight  medical  treatment  ability  and  flight  speed.  Constraints  (16)  assign  one  type  A  air  asset  in  each  of  the  first  1, . . . ,  sa 
positions  in  the  preference  list  to  casualties  of  type  Ca  or  Cab  ■  Casualties  of  type  Cba  are  assigned  all  of  the  type  B  air 
assets  first  in  the  preference  list  before  the  type  A  air  assets  are  assigned  in  the  preference  list  as  seen  in  constraints  (17). 
Likewise,  only  type  B  air  assets  respond  to  Cb  casualties  and  Cba  casualties  prefer  type  B  air  assets.  Constraints  (18) 
assign  one  type  B  air  asset  in  each  of  the  first  1, . . . ,  sb  positions  of  the  preference  list  to  casualties  of  type  Cb  or  Cba- 
Casualties  of  type  Cab  are  assigned  all  of  the  type  A  air  assets  in  the  preference  list  first  before  the  type  B  air  assets  are 
assigned  in  the  preference  list  as  seen  in  constraints  (19).  Constraints  (20)  -  (23)  assign  each  air  asset  in  the  system  to 
only  one  position  in  the  preference  list.  In  other  words,  a  specific  air  asset  cannot  be  both  the  first  and  third  preferred 
air  asset  for  the  same  casualty  location.  Constraints  (24),  (25),  and  (26)  enforce  minimum  pre-specified  thresholds  9C 
of  the  expected  coverage  for  the  non-objective  casualty  types  c  €  Cab  U  CbaC  Cb-  The  minimum  levels  of  service  are 
typically  defined  by  a  division  surgeon  cell  or  senior  medical  planners.  Constraints  (27)  and  (28)  locate  sa  type  A  and  Sb 
type  B  air  assets  at  MTFs  in  the  system,  respectively.  Constraints  (29)  -  (32)  make  sure  that  only  air  assets  located  at 
open  MTF  facility  locations  are  utilized  in  construction  of  the  preference  list.  Constraints  (33)  compute  and  balance  the 
type  A  workloads  within  a  pre-specified  < 5A  of  the  average  workload.  In  a  similar  manner,  constraints  (34)  compute  and 
balance  the  type  B  workloads.  Constraints  (35)  and  (36)  enforce  contiguity  amongst  each  type  A  and  type  B  air  asset’s 
set  of  first  priority  casualty  locations.  Two  locations  are  contiguous  to  each  other  if  the  two  locations  differ  by  one  unit 
in  either  the  x  or  y  direction.  Lastly,  constraints  (37)  -  (38)  and  (39)  -  (40)  require  the  preference  list  and  the  locating 
variables  to  be  binary,  respectively.  In  the  next  section,  we  present  an  example  to  apply  the  BLP  model  to  a  military 
medical  system. 


5  Computational  Example 

We  now  present  a  computational  example  to  illustrate  how  to  locate  two  classes  of  air  assets  and  assign  them  to  multiple 
types  of  casualties  in  a  dispatch  preference  list.  In  the  ongoing  stability  operations  in  support  of  Operation  Enduring 
Freedom  (OEF)  in  Afghanistan,  coalition  forces  from  the  North  Atlantic  Treaty  Organization  (NATO)  leverage  both  U.S. 
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Army  and  U.S.  Air  Force  helicopters  in  support  of  combat  troops  conducting  operations  on  the  ground.  Specifically,  the 
U.S.  Army  air  assets  (UH-60A/L  and  HH-60  MEDEVAC)  and  the  U.S.  Air  Force  air  asset  (HH-60G  Pave  Hawk)  are 
primarily  responsible  for  the  aeromedical  evacuation  of  casualties  [2] .  The  HH-60  MEDEVAC  platform  -  referred  to  as  air 
asset  type  A  has  been  specifically  modified  with  a  state-of-the-art  medical  interior  to  accommodate  acute  care  casualties 
and  has  a  cruising  speed  of  278  km/hour.  The  HH-60G  Pave  Hawk  air  asset  -  referred  to  as  air  asset  type  B  is  a  slower 
air  asset  with  a  cruising  speed  of  240  km/hour  and  does  not  have  the  in-flight  medical  capabilities  to  accommodate 
acute  care  casualties  [22],  In  this  paper,  we  assume  that  both  type  A  and  type  B  air  assets  fly  single-ship  (i.e. ,  solo) 
during  MEDEVAC  operations,  in  contrast  to  tandem  operations  in  which  type  A  air  assets  must  be  accompanied  by  a 
security/escort  asset  due  to  risk  of  enemy  action.  Knowing  where  to  locate  both  type  A  and  type  B  air  assets  and  how 
to  dispatch  them  to  battle  triaged  casualties  will  help  increase  survivability  of  the  most  urgent  casualties. 

We  apply  our  air  MEDEVAC  asset  optimization  model  to  NATO  Regional  Command  -  South  (RC-South),  one  of 
the  four  regional  commands  which  span  several  Afghanistan  provinces,  to  illustrate  how  air  asset  management  can  lead 
to  more  effective  military  medical  systems.  Figure  1  depicts  the  area  of  RC-South  divided  into  seventy-five  distinct 
casualty  locations.  Twenty- five  potential  sites  for  medical  treatment  facilities  (MTF)  are  distributed  throughout  RC- 
South,  represented  by  the  red  crosses  in  Figure  1.  The  base  case  example  solved  in  this  section  will  optimally  locate  ten 
type  A  air  assets  and  four  type  B  air  assets  at  MTFs  within  RC-South.  Note  that  the  casualty  data  (from  2001  through 
2012)  was  one-time  extracted  from  the  Defense  Casualty  Analysis  System  (DCAS). 


Fig.  1:  25  Potential  MTF  sites  in  NATO  Regional  Command-South  (represented  by  red  crosses) 

The  following  organization  of  casualty  types  is  used  to  be  consistent  with  the  notation  of  sections  3  and  4.  Casualty 
type  Ca  contains  the  CAT  A  casualties,  Cab  contains  the  CAT  B  casualties,  and  Cb  contains  the  CAT  C  casualties. 
Note  that  there  are  no  Cba  call  types  in  this  example. 

We  next  derive  the  BLP  model  inputs.  The  type  A  service  times  t£c  (hours)  are  computed  by  dividing  the  Euclidean 
distance  (kilometers)  between  a  MTF  location  i  G  I  a  and  a  casualty  location  j  £  J  by  the  type  A  air  asset  speed 
(kilometers/hour).  The  type  B  service  times  t,^c  are  computed  similarly  except  that  the  service  time  calculation  is 
multiplied  by  two  since  medical  treatment  does  not  begin  until  a  casualty  returns  back  to  the  MTF  (in  contrast  to  the 
medical  treatment  provided  at  the  casualty  location  with  a  type  A  air  asset).  The  distribution  of  casualty  type  is  as  follows: 
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CAT  A  75%,  CAT  B  20%,  and  CAT  C  5%.  Further,  the  total  casualty  arrival  rate  A  =  3.0  and  each  casualty  location  j  £  J 
has  equal  chance  of  a  casualty  event.  Evacuation  of  a  CAT  A  casualty  by  a  type  A  air  asset  results  in  an  expected  coverage 
between  0.8  1.0,  depending  on  how  many  minutes  have  elapsed  up  to  the  threshold  of  one  hour.  Therefore,  we  use  an 

exponential  decreasing  curve  to  model  the  rapid  deteriorating  health  of  CAT  A  casualties  after  one  hour.  Similar  translated 
exponential  decreasing  curves  are  used  for  the  evacuation  of  CAT  B  and  CAT  C  casualties  to  reflect  time-standards  of  four 
hours  and  twenty-four  hours,  respectively.  We  use  a  time-standard  of  six  hours  for  CAT  C  casualties  to  encourage  the  rapid 
evacuation  of  all  casualty  types.  Rf-C  and  RBC  depend  on  the  expected  type  A  and  type  B  air  asset  service  time  between 
casualty  location  j  and  air  asset  location  i  (denoted  t  in  the  equations  below)  as  well  as  the  casualty  type  c  £  C.  R.fjC  and 
RBC  and  are  given  by  the  following  equations.  For  c  £  CAT  A:i?ijc  =  —  0.2£  +  l  if  t  <  1  and  e_1-6(t_1)  otherwise  ,  For  c  £ 
CAT  B :RijC  =  —  0.05t+l  if  t  <  4  and  e-1'6^-4^  otherwise  ,  and  For  c  £  CAT  C :i?,yc  =  —  0.05t  +  l  if  t  <  6  and  e-1-6^-6) 
otherwise. 

We  perform  all  computations  on  laptop  with  a  AMD  A6  2.10  GHz  processor  and  8GB  RAM,  leveraging  Python  2.6 
for  parameter  calculations  and  Gurobi  5.6.2  for  optimization  of  the  BLP  model.  The  following  algorithm  is  executed  to 
solve  the  optimization  model.  Note  that  each  optimization  iteration  contains  an  instance  of  the  approximate  Hypercube 
algorithm.  Optimize  the  BLP  using  the  most  recent  values  of  the  following  system  statistics:  1)  the  average  busy  proba¬ 
bilities  rA  and  rB ,  2)  the  independence  correction  factors  Qa  and  Qb,  and  3)  the  offered  load  imbalance  thresholds  5a 
and  6b-  Optimize  a  new  instance  of  the  problem  until  the  workload  of  each  type  A  and  type  B  air  assets  is  less  than  5% 
more  than  the  average  air  asset  workload.  After  each  optimization  iteration,  the  current  solution  contains  the  location 
of  the  air  assets  and  the  dispatch  preference  list.  Using  the  dispatch  preference  list,  execute  the  approximate  Hypercube 
algorithm  until  the  max  change  in  type  A  and  type  B  busy  probabilities  is  less  than  0.001.  The  system  performance 
statistics  ( rA ,  rB ,  Qa(sa ,  pA ,  k),  Qb(sb,  pB ,k))  are  updated  and  passed  as  input  to  the  next  instance  of  the  BLP  model. 
The  offered  load  imbalance  thresholds  5a  and  5b  are  tightened  to  the  new  standard  deviation  of  the  type  A  and  type  B 
busy  probabilities.  For  the  base  case  example,  the  optimization  algorithm  converged  in  four  iterations.  The  approximate 
Hypercube  algorithm  converged  in  six  or  fewer  iterations  within  each  iteration  of  the  optimization  algorithm.  The  base 
case  example  corresponds  to  a  BLP  with  67,757  constraints  and  52,602  variables.  Solving  the  optimization  algorithm  on 
the  base  case  ten  times  resulted  in  the  following  run  time  metrics:  a  minimum  CPU  time  of  908  seconds  (15.1  minutes), 
a  maximum  CPU  time  of  1,132  seconds  (18.9  minutes),  and  an  average  CPU  time  of  1,001  seconds  (16.7  minutes). 

Table  1  presents  the  ten  type  A  and  four  type  B  air  asset  locations  selected  by  the  BLP  model  as  well  as  the 
corresponding  busy  probabilities.  MTF43  is  the  only  location  with  both  a  type  A  and  a  type  B  air  asset.  The  type  A 
busy  probabilities  range  from  the  low  of  0.051  (MTF07)  to  the  high  of  0.103  (MTF25).  The  type  B  busy  probabilities 
range  from  the  low  of  0.224  (MTF37)  to  the  high  of  0.358  (MTF43).  All  air  asset  busy  probabilities  are  balanced  to  be 
no  greater  than  5%  more  than  the  respective  type  A  or  type  B  average  busy  probability. 

Table  2  presents  the  independence  correction  factors  for  type  A  and  type  B  air  assets.  The  type  A  factors  quantify 
the  correction  to  the  probability  of  obtaining  k  busy  type  A  air  assets  followed  by  an  available  type  A  air  asset  when 
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Table  1:  Location  and  Busy  Probabilities  for  Type  A  and  Type  B  Air  Assets  in  the  Base  Case  Example 


The  location  of  the  sa  —  10  type  A  air  assets 

MTF25 

MTF22 

MTF28 

MTF07 

MTF46  MTF61 

MTF64 

MTF43 

MTF10 

MTF58 

Type  A  air  asset  busy  probabilities  rA 

0.103 

0.058 

0.098 

0.051 

0.068  0.086 

0.076 

0.061 

0.081 

0.070 

Average  type  A  air  asset  busy  probability  rA 

0.075 

The  location  of  the  sb  —  4  type  B  air  assets 

MTF43 

MTF73 

MTF37 

MTF  19 

Type  B  air  asset  busy  probabilities  rf 

0.358 

0.339 

0.224 

0.320 

Average  type  B  air  asset  busy  probability  r B 

0.310 

assuming  that  air  assets  operate  independently.  Note  that  Qa(sa,PA,  0)  =  Qb(sb,  pB ,  0)  =  1.0  because  the  probability 
of  obtaining  an  available  air  asset  when  no  air  assets  are  busy  is  equal  to  one. 


Table  2:  Independence  Correction  Factors  for  Type  A  and  Type  B  Air  Assets  in  the  Base  Case  Example 


QA{sA,PA,k ) 

k  =  0 

1.0 

k=  1 
0.907 

k  =  2 
0.878 

k =3  k= 4  k= 5  k =6  k  = 7  k = 8  k= 9 

0.907  0.995  1.155  1.411  1.805  2.405  3.318 

Qb{sb,PB  ,k) 

k  =  0 

1.0 

k  =  1 
0.856 

k  =  2 
0.807 

k  =  3 

0.830 

The  objective  function  of  the  BLP  is  the  expected  fraction  of  calls  covered  when  dispatching  type  A  air  assets  to 
CAT  A  casualties.  It  is  desirable  to  know  the  performances  changes  in  the  system  as  additional  type  A  air  assets  are 
introduced  or  removed  from  the  system.  Figure  2  presents  a  sensitivity  analysis  on  the  expected  coverage  of  CAT  A 
casualties  with  respect  to  the  number  of  type  A  air  assets  in  the  system  as  well  as  the  call  arrival  rate  A.  Figure  2  can  be 
used  by  decision-makers  to  identify  how  many  type  A  air  assets  are  needed  for  a  given  benchmark.  For  example,  if  the 
desired  system  performance  is  90%  expected  coverage,  then  six  type  A  air  assets  will  be  enough  when  the  call  volume  A  is 
between  2.5  and  3.5  calls  per  hour.  However,  to  achieve  95%i  expected  coverage  in  the  same  call  volume  range,  nine  type 
A  air  assets  are  required.  In  figure  2,  there  is  a  diminishing  return  on  the  objective  function  value  with  respect  to  the 
number  of  type  A  air  assets.  For  example,  when  A  =  3.0,  consider  the  large  increase  in  expected  coverage  when  sa  =  2 
(53%)  and  when  =  5  (89%)  to  the  small  increase  when  s,4  =  8  (94%)  and  when  sa  =  10  (95%). 

Figures  3(a)  -  3(c)  illustrate  the  contiguity  of  the  1st  priority  response  districts  created  for  CAT  A,  CAT  B,  and  CAT 
C  casualties,  respectively.  Figures  3(a)  and  3(b)  show  type  A  air  assets  because  the  first  priority  for  CAT  A  and  CAT 
B  casualties  must  be  a  type  A  air  asset.  The  BLP  model  identifies  the  set  of  type  A  air  asset  locations  that  maximize 
the  fraction  of  CAT  A  casualties  covered  within  60  minutes.  The  set  of  locations  selected  must  also  balance  air  asset 
workloads  while  considering  busy  probabilities  and  accounting  for  air  assets  that  are  dependent.  Likewise,  figure  3(c) 
shows  type  B  air  assets  because  the  first  priority  for  each  CAT  C  casualty  must  be  a  type  B  air  asset. 

The  objective  function  used  in  this  computational  example  encourages  dispatch  of  the  closest  air  asset  to  a  casualty 
incident.  Dispatching  the  closest  air  asset  leads  to  the  shortest  service  time  and  subsequently  the  largest  expected  coverage. 
However,  due  to  the  unavailability  of  air  assets  in  a  crowded  system,  it  is  sometimes  optimal  to  ration  the  closest  available 
air  asset.  Table  3  presents  the  difference  in  objective  function  value  between  optimality  and  the  heuristic  policy  of  dispatch 
the  closest  air  asset  to  a  casualty.  To  implement  the  closest  air  asset  heuristic,  sa  type  A  and  sb  type  B  air  assets  must 
be  selected  for  use  in  the  priority  list.  Therefore,  a  rank  ordering  of  all  potential  MTF  sites  is  created  based  upon  the 
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Fig.  2:  Sensitivity  Analysis  of  the  Objective  Function  for  the  Number  of  Type  A  Air  Assets  and  Call  Arrival  Rate  A. 


Source  for  ISAF  Regions:  http://www.nato.int/isaf 


(a)  1st  Priority  Response  Districts  of  Type  A  Air  Assets  (b)  1st  Priority  Response  Districts  of  Type  A  Air  Assets 
Assigned  to  CAT  A  Casualties  Assigned  to  CAT  B  Casualties 


Source  for  ISAF  Regions:  http://www.nato.int/isaf 


(c)  1st  Priority  Response  Districts  of  Type  B  Air  Assets 
Assigned  to  CAT  C  Casualties 

Fig.  3:  Contiguous  1st  priority  dispatching  assignments  for  type  A  and  type  B  air  assets  in  the  base  case  example.  Casualty 
locations  that  comprise  an  air  asset’s  1st  priority  district  are  the  same  color. 


average  distance  to  all  seventy  five  casualty  locations.  Once  the  air  asset  locations  are  known,  the  priority  list  is  created 
based  on  proximity  between  each  casualty  location  and  each  asset  location.  To  evaluate  the  objective  function  value  of 
the  heuristic,  to  locating  and  dispatching  variables  are  set  to  1.0,  the  system  input  parameters  are  calculated,  and  the 
instance  is  solved  in  Gurobi.  Table  3  shows  that  the  heuristic  of  dispatching  the  closest  air  asset  performs  very  well.  The 
largest  difference  in  objective  function  value  between  optimality  and  the  heuristic  over  the  range  of  and  sb  values 
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considered  is  3.99%.  In  some  scenarios,  such  as  sa  =  2  and  sb  =  4,  the  heuristic  policy  is  optimal.  Note  that  in  Table  3 
the  base  case  value  of  A  =  3  is  used. 

Table  3:  Optimal  Solution  Versus  Closest  Server  Heuristic  Solution  with  Additional  Type  A  Air  Assets 


sa 

SB 

Heuristic  Obj  Value 

Optimal  Obj  Value 

Obj  Improvement  (percent) 

2 

4 

0.533 

0.533 

0% 

3 

4 

0.710 

0.710 

0% 

4 

4 

0.805 

0.834 

3.60% 

5 

4 

0.853 

0.887 

3.99% 

6 

4 

0.884 

0.912 

3.17% 

7 

4 

0.906 

0.926 

2.21% 

8 

4 

0.916 

0.937 

2.29  % 

9 

4 

0.922 

0.944 

2.39% 

10 

4 

0.933 

0.949 

1.71% 

The  comparison  of  optimality  to  the  dispatch  closest  heuristic  is  examined  further  in  Table  4.  The  call  arrival  rate 
is  varied  from  A  =  1  to  10.  Under  this  range  of  values,  the  heuristic  of  dispatching  the  closest  air  asset  again  performs 
very  well.  The  largest  difference  in  heuristic  objective  function  value  from  optimality  is  1.71%.  Note  that  in  Table  4  the 
base  case  values  of  sa  =  10  and  sb  =  4  are  used.  These  improvements,  while  small  in  absolute  value,  are  large  relative 
to  other  methodological  improvements  in  this  domain. 


Table  4:  Optimal  Solution  Versus  Closest  Server  Heuristic  Solution  with  Varying  Call  Arrival  Rate  A 


Call  Arrival  Rate  A 

Heuristic  Obj  Value 

Optimal  Obj  Value 

Obj  Improvement  (percent) 

1 

0.941 

0.957 

1.70% 

3 

0.933 

0.949 

1.71% 

6 

0.918 

0.932 

1.53% 

10 

0.908 

0.911 

0.33% 

This  computational  example  illustrates  how  to  locate  two  classes  of  air  assets  and  construct  an  air  MEDEVAC  asset 
dispatch  preference  list  for  all  casualty  locations.  The  base  case  example  developed  in  this  section  illustrates  how  to  locate 
and  dispatch  10  type  A  and  4  type  B  air  assets  to  maximize  the  expected  coverage  of  CAT  A  casualties  within  sixty 
minutes,  while  balancing  the  air  asset  workload.  A  dispatch  preference  list  is  valuable  to  senior  military  decision-makers 
and  medical  planners  to  identify  preferred  back-up  air  assets  in  a  crowded  military  MEDEVAC  system. 


6  Concluding  Remarks 

In  this  paper,  we  introduced  a  binary  linear  programming  model  that  optimally  locates  two  classes  of  military  air 
MEDEVAC  assets  and  assigns  air  assets  to  casualties  in  a  dispatch  preference  list.  Queuing  system  performance  statistics, 
including  busy  probabilities  and  independence  correction  factors,  are  derived  for  air  assets  and  are  leveraged  as  inputs  to 
the  BLP  model.  Accounting  for  the  queuing  dynamics  in  military  medical  systems  results  in  a  more  realistic  optimization 
model  in  which  air  assets  are  not  always  available  for  service.  The  problem  of  military  air  MEDEVAC  asset  management 


Optimally  Locating  and  Dispatching  Two  Classes  of  Military  MEDEVAC  Air  Assets 


17 


is  important  because  military  medical  systems  save  lives  by  responding  to  multiple  types  of  casualty  incidents.  By 
providing  prompt  en  route  medical  care  to  the  most  urgent  casualties,  the  likelihood  of  survival  is  increased.  The  effective 
and  efficient  medical  evacuation  of  casualties  serves  to  keep  troop  morale  high  while  maintaining  a  healthy  and  sustained 
fighting  force. 

The  model  in  this  paper  can  be  extended  to  consider  more  than  two  types  of  air  assets.  The  resulting  spatial  queuing 
model  and  BLP  can  be  derived  using  methods  that  are  analogous  to  those  in  this  paper. 

This  paper  also  can  be  extended  to  account  for  batch  arrivals  of  casualties  to  the  system  where  the  volume  of 
casualties  requires  the  dispatch  of  more  than  one  air  asset.  Also,  the  initial  triage  of  casualties  is  not  perfect  and  can 
lead  to  misclassification  of  the  severity  of  a  casualty  incident.  The  logistics  of  military  medical  systems  are  even  more 
complicated  when  there  are  multiple  air  assets  and  then  the  accuracy  of  triage  affects  dispatch  decisions.  Also,  in  some 
cases,  air  assets  must  be  assigned  to  and  wait  for  security /escort  assets.  An  optimization  model  to  identify  which  air 
asset  and  security /escort  asset  pair  to  dispatch  to  casualty  incidents  to  minimize  MEDEVAC  response  time  of  urgent 
casualties  would  be  of  value  to  military  decision-makers.  Work  is  currently  in  progress  to  address  these  issues. 
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